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Abstract 


Free  viscously  damped  vibrations  of  linear  discrete  structural  systems  are 
studied.  The  amount  of  damping  varies  among  the  various  structural  elements 
of  the  system  resulting  in  several  critical  damping  possibilities.  A  general 
method  is  developed  far  determining  the  "critical  damping  surfaces"  of  a 
system.  These  surfaces  represent  the  loci  of  combinations  of  damping  values 
corresponding  to  critically  damped  motions,  and  thus  separate  regions  of 
partial  or  complete  underdamping  from  those  of  overdamping.  The  dimension  of 
a  critical  damping  surface  is  equal  to  the  number  of  independent  amounts  of 
damping  present  in  the  system.  The  determination  of  the  surface  point  cor¬ 
responding  to  equal  amounts  of  damping  is  considerably  simplified  for  systems 
which,  on  the  assumption  that  all  amounts  of  damping  are  equal,  possess  a 
damping  matrix  of  the  Rayleigh  type.  Three  examples  presented  in  detail  illus¬ 
trate  the  proposed  technique  and  some  of  the  important  characteristics  of 


critical  damping  surfaces. 
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Introduction 

The  amount  of  damping  in  a  linear  structural  dynamic  system  is  usually 
expressed  as  a  percentage  of  the  critical  or  the  modal  critical  damping  (e.g. 
[1]).  It  is  then  easy  to  determine  if  the  structure  is  underdamped  or  over¬ 
damped  and  consequently  to  control  the  response  by  appropriately  varying 
the  amount  of  structural  damping.  However,  this  procedure  is  possible  only 
when  the  damping  is  the  same  everywhere  in  the  structure  or  is  given  in  modal 
form,  in  which  cases  the  problem  of  determining  the  amount  of  critical  damp¬ 
ing  is  essentially  one-dimensional  and  presents  no  difficulties. 

However,  for  a  linear  structural  dynamic  system  with  viscous  damping 
varying  among  its  elements  the  problem  of  determining  critical  damping  becomes 
much  more  difficult  since  there  are  many  critical  and  partially  critical  damp¬ 
ing  possibilities  corresponding  to  different  critical  element  damping  value 
combinations.  The  importance  of  this  problem  lies  in  the  fact  that  introduces 
the  possibility  of  controlling  the  dynamic  response  more  easily  and  with  a 
greater  flexibility  by  differently  varying  the  damping  of  a  number  of  elements. 

Critical  damping  analyses  are  usually  performed  in  connection  with  free 
vibration  studies  which  aim  at  determining  natural  frequencies  and  modal  shapes. 
Necessary  and  sufficient  conditions  under  which  discrete  damped  linear  dynamic 
systems  possess  classical  normal  modes  have  been  established  by  Caughey  [2] 
and  Caughey  and  O' Kelly  [3], 

For  a  general  discrete  viscously  damped  dynamic  system  with  n  degrees  of 


freedom,  the  free  vibration  problem  primarily  consists  of  determining  the 
jrogts  of  ^he^determinantal  or  characteristic  equation  of  the  system.  These 
iflcarl  b^  iWl  V  ex  or  purely  imaginary  and  thus  characterize  the  motion  as 
r-  ,  overdamped,  unlerdamped  or  undamped,  respectively.  This  is  what  is  usually 
mentioned  in  the  literature  (e.g.  [1])  without  consideration  of  the  case  of 
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crltical  damping  or  the  case  of  coexistence  of  all  kinds  of  the  above  roots. 

The  only  exception  appears  to  be  Crafton  [4]  who  considers  two  degrees  of 
freedom  spring-mass  systems  with  two  different  amounts  of  viscous  damping, 
and  determines  the  (in  general)  complex  roots  of  the  characteristic  equation, 
omitting,  however,  a  critical  damping  analysis.  The  complex  roots  are  ob¬ 
tained  by  rather  cumbersome  conformal  mapping  techniques  in  [4],  but  even 
If  more  efficient  numerical  techniques  are  employed  (e.g.  [5]),  this  method 
of  characterization  of  the  motion  is  not  convenient  for  design  purposes, 
since  It  requires  a  complete  root  determination  for  every  combination  of 
element  damping  values  assumed.  Bishop  and  Johnson  [6],  on  the  basis  of  a 
particular  example  of  a  two  degrees  of  freedom  system  with  only  one  viscous 
damper,  studied  the  behavior  of  the  frequencies  of  the  system  as  damping  in¬ 
creased  form  zero  to  Infinity. 

The  present  paper  presents  a  detailed  study  of  the  effect  of  damping 
distributed  in  an  arbitrary  manner  throughout  the  structural  elements  on  the 
free  motion  of  the  structure.  A  general  method  is  first  proposed  to  determine 
the  "critical  damping  surfaces"  of  a  viscously  damped  linear  discrete  dynamic 
system;  these  are  the  loci,  in  "damping  space,"  of  amounts  of  damping  leading 
to  critically  damped  motions.  Three  examples  are  then  presented  in  detail 
which  illustrate  the  proposed  method  and  lead  to  some  interesting  and  unexpected 


results. 
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Free  Vibrations  of  Damped  Oiscrete  Systems 

Consider  a  viscously  damped  linear  discrete  system  of  n  degrees  of  free¬ 
dom  In  free  motion  governed  by  the  matrix  equation 

[M]{x}  +  [C](x>  +  [K]{x}  3  {0},  (1) 

where  [M],  [C]  and  [K]  are  symmetric  inertia,  damping  and  stiffness  n  x  n 
matrices,  respectively,  {x}  3  (x(t)}  is  the  n  x  1  displacement  vector  arid 
dots  indicate  differentiation  with  respect  to  time  t. 

Assuming  a  solution  of  the  form 

{x>  *  ($}eXt  ,  (2) 

one  reduces  (1)  to  the  eigenvalue  problem 

[A2[M]  +  X[C]  +  [*]](*}  3  [D( A) ]{$}  3  {0}  ,  (3) 

which  has  nontrivial  solutions  if,  and  only  if, 

det[D(A)]  =  | 0(A) |=0.  (4) 

Equation  (4)  is  the  characteristic  or  determinantal  equation,  an  algebraic 
equation  of  order  2n  In  A  .  The  roots  of  (4)  can  be  real ,  imaginary  or  com¬ 
plex  and  occur  in  pairs  of  the  general  form 

A  «  -u±4»  ,  (5) 

where  \i  Is  nonnegative  real,  with  its  negative  sign  indicating  stability 
of  the  system,  and  if/  Is  a  nonnegative  real  or  imaginary  number.  A  real  root  A, 
which  should  be  negative  for  stability,  indicates  either  overdamping  (for 
which  ijKy),  or  critical  damping  (for  which  ip=0) ;  both  cases  correspond  to 
aperiodic  decaying  motions.  A  complex  root  A,  for  which  is  an  imaginary 
number  of  the  form  jw  ( ja/-T) ,  Indicates  underdamping  corresponding  to  an 
oscillatory  decaying  motion  with  a  damped  natural  frequency  w.  For  an  un¬ 
damped  system  for  which  [C]  3  [0],  all  the  roots  of  (4)  are  given  by  (5)  with 
y=0  and  <J>=jw0  and  correspond  to  periodic  motions  of  natural  frequencies  wQ  . 

If  some  (all)  of  the  roots  of  (4)  correspond  to  overdamping,  underdamping  or 
critical  damping,  the  system  is  called  partially  (completely)  overdamped. 


underdamped  or  critically  damped,  respectively. 

Once  the  roots  (£  *  l,2,...,2n)  of  (4)  have  been  determined,  one  can 
solve  (3)  for  the  modal  shapes  which  are,  in  general,  complex.  Thus, 

the  solution  of  (1)  takes  the  general  form 

{x}  =  [$]  {QeXt}  (6) 

where  [$]  is  the  n  x  2n  modal  matrix  and  Q  are  arbitrary  real  or  complex 
constants  to  be  determined  from  the  initial  conditions  of  the  problem.  A 
necessary  and  sufficient  condition  for  a  passive  viscously  damped  linear 
discrete  system  to  possess  classical  normal  modes,  i.e.,  a  complete  set  of 
real  orthonormal  eigenvectors,  is  that  [3] 

[C][M]-l[K3  =  ,  (7) 

In  which  case  these  classical  normal  modes  are  identical  with  the  normal 
modes  of  the  undamped  system. 
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Critical  Oampinq  Surfaces 

Consider  a  n  degrees  of  freedom  linear  damped  dynamic  system  character¬ 
ized  by  the  matrices  [M],  [K],  and  [C],  where  the  elements  of  [C]  are  combina¬ 
tions  of  the  m  different  damping  coefficients  ck  (k  =  l,2,...,m)  of  the  various 
elements  of  the  system.  The  number  m  can  be  greater,  equal  or  smaller  than 
n.  For  partial  (complete)  overdamping  or  critical  damping  some  (all)  of  the 
roots  of  (4)  are  of  the  form 

X  -  -b  ,  b  >  0  (8) 

and  in  that  case  (4)  becomes 


|b2[M]  -b[C]  +  [K]|  =  |D(b,ck)|  =0.  (9) 

In  the  m-dimensional  space  with  coordinates  cfc  (k  =  l,2,...,m),  equation  (9) 
represents  a  family  of  m-dimensional  surfaces  S^.  Each  member  of  the  family 
corresponds  to  overdamping  or  critical  damping  and  is  characterized  by  its  own 
value  of  b,  which  is  a  function  of  the  ck's.  The  problem  is  to  determine  that 
b  which  corresponds  to  the  "critical  damping  surface,"  i.e.,  the  locus  of 
combinations  of  ck  leading  to  critically  damped  motion  and  thus  separating 
regions  of  partial  or  complete  underdamping  from  those  of  overdamping.  There 
are  actually  q  (q  *  n)  critical  damping  surfaces  since  there  are  at  most  as 
many  partial  critical  damping  possibilities  as  the  number  of  the  pairs  of  roots 
X  in  (5)  with  i|i  *  0,  Since  critical  damping  represents  the  threshold  between 
overdamping  and  underdamping,  one  can  conclude  that  among  the  Sb  surfaces,  the 
critical  surface  Scr  is  the  one  for  which  the  damping  will  be  a  minimum,  i.e.. 


(d/db)  |D(b,ck) |*(3/3b) |D(b,ck) |+E(3/ack)  |0(b,ck)  | (3ck/3b) 
(3ck/3b)  =0  ,  k  *  1 ,2,.. .,m. 


0 


J 


(10) 


An  alternative  derivation  of  (9)  and  (10)  is  the  following:  Consider 
equation  (4)  with  X  *  -R+jl,  (R>0,  I*real,  j-/-l),  i.e.. 


A(X)  *  A(-R,jl)  -  |D(-R+jI)|  *  0  . 


(ID 
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Expansion  of  A(A)  in  Taylor  series  about  the  point  AQ  *  ( -R,0)  yields 

A(-R,jl)  *  A(-R,0)  +  (jl)  £  l(.Rj0)  +  jM'^fx^-R.O)  + 

+  Tl^1  !(-R,0)  +  •  ••  '  (12* 

which  in  view  of  (11)  reduces  to 

a("R»o)  -  2!  12  Ip-  I(-r#o)  +  *  0  ’ 

,  03) 

~  I  _  I  t2  9A,  |  +  « 

3X  *(-R,0)  3!  1  aF  l(-R,0)  •**  0  *  J 

For  critical  damping,  for  which  1=0,  equations  (13)  lead  to  equations 

A(-R,0)  -  0  , 

04) 

(54/3A)I(-R,Q)  ■  °>  J 

which  are  the  same  as  (9)  and  (10). 

In  principle,  one  can  solve  (10)  for  b  and  obtain  its  critical  value  bcp 
as  a  function  of  the  c^’s.  Thus,  the  equation  of  critical  damping  surfaces 
will  be  given  by  (9)  with  b  *  b_„,  i.e.,  by 

Ci 

|b*r[M]  -  bcr[C]  +  [K]|  *  0  .  (15) 

In  practice,  however,  the  nonlinear  system  of  (9)  and  (10)  has  to  be  solved 
numerically.  Thus,  provided  that  differentiations  in  (10)  can  be  done  ana¬ 
lytically,  one  determines  numerically  a  finite  number  of  groups  of  m+1  values 
for  bcr  and  c^  satisfying  (9)  and  (10)  simultaneously,  and  one  is  then  able 
to  represent  critical  damping  surfaces  in  the  m-dimensional  c^  space  as  sets 
of  points  with  coordinates  c^.  Eventhough  the  above  method  of  determining 
critical  damping  surfaces  Is  quite  general  and  applicable,  in  principle,  to  a 
n  degrees  of  freedom  system,  the  fact  that,  to  the  authors  knowledge,  there  is 
no  presently  available  efficient  numerical  treatment  of  the  differentiations  in 
(10),  limits  the  range  of  applicability  of  the  method  to  small  order  systems. 
However,  the  particular  point  ct  *  c2  *  ...  *  cm  *  c  of  a  critical  damping 


surface  can  be  very  easily  obtained  for  the  special  case  in  which  [CL 

vc 

is  of  the  Rayleigh  type,  i.e.,  of  the  form 

[C]c  *  a  [M]  +  a  [K]  ,  (16) 

ck  c  i  2 

where  a  and  a  are  constants.  For  this  case  (3)  becomes 
1  2 

[(X2+a  A) CM3  +  (1+a  A)[K]]{$}  =  {0},  (17) 

x  2 

indicating  that  there  exist  classical  normal  modes  and  that  for  the  i ^ 

(i  *  l,2,...,n)  mode  with  underdamping 

(X2  +  aiXi)/(l  +  a^)  »  -  w^  .  (18) 

Equation  (18)  can  be  solved  for  X.  and  give 

1* 


X.  3  -b  .±jw .  , 

i  t  t 


where 


b.  -  l/2(a  +a  w2.)  ,  w?  =  w2.  -  b2  , 

1  12  01  l  0^  t 

,  --  '  (20) 

*.  *  *,  "Si  ‘  0  -  J  ^  ■  J 

It  Is  apparent  then,  in  view  of  (19)  and  (20),  that  at  critical  damping  one 

has  w.  *  0,  or 

"icr  '  "oi  •  <2” 

with  a^  and  a^  satisfying  the  relation 

Vr  “Si  -  2"oi  +  aicr  ■  0  •  <22> 

where  a  and.  a  are  their  values  at  critical  damping.  Equations  (21)  and  (22) 
xcr  2cr 

therefore  represent  the  solution  of  (10)  for  the  particular  surface  point 
c^  *  c2  a  3  cm*c  and  for  systems  obeying  (16).  This  solution  Is  unique  if 

one  of  the  a  and  a  is  zero.  Once  the  b. 's  have  been  computed  from  (21), 

a  numerical  evaluation  of  the  left  hand  side  of  (15)  for  a  sequence  of  values 
of  c  leads  to  the  determination  of  the  common  value  c  of  the  c^'s  which  satisfies 
(15)  for  every  i.  One  important  result  of  this  analysis  Is  that,  as  equation 
(21)  clearly  demonstrates,  for  systems  obeying  (16),  one  has  q»n,  i.e.,  as  many 
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critical  damping  surfaces  as  there  are  degrees  of  freedom  in  the  system. 

Three  illustrative  examples  follow. 

Example  1 

For  the  familiar  one  degree  of  freedom  spring-dashpot-mass  system  the 
equation  of  free  motion  is 

x  +  2gx  +  w2  *  0  ,  (23) 

where  the  measure  of  the  amount  of  damping  8  and  the  undamped  natural  frequency 
wQ  are  given  in  terms  of  the  mass  m,  the  stiffness  K  and  the  damping  coefficient 
c  of  the  system  as 

8  *  c/2m  ,  wQ  =  A/m  .  (24) 

Here  m=nsl  so  that  there  exists  only  one  critical  damping  surface,  of  dimension 

one,  i.e. ,  a  point.  The  condition  (9)  for  overdamping  or  critical  damping  is 

b2-  2bg  +  w2  3  0  .  (25) 

Use  of  (10)  provides  b  =  8  which,  in  conjunction  with  (25),  leads  to 

Ci  c  * 

8cr  *  wq  ,  (26) 

which  is  the  well  known  condition  of  critical  damping  for  this  system.  Notice 

that  (26)  could  also  have  been  obtained  from  (20),  (21)  and  (24)  by  realizing 

that  c  ■  (c/m)m  =  a  m  ,  in  which  case  b__  *  (l/2)a  =  c  /2m  =  8__  *  w  . 

i  Ci  iCi  c*  c«  o 

Example  2 

Consider  the  two  degrees  of  freedom  system  of  Fig.  1  consisting  of  a  vari¬ 
able  torsional  stiffness  shaft  with  two  flywheels  and  one  viscous  damper  and 
being  in  free  torsional  vibration.  This  system  has  been  taken  from  reference 
[6]  (p.  507)  which  studies  the  way  in  which  the  free  motion  of  the  system  changes 
as  the  damping  is  increased  from  zero  to  infinity  by  determining  the  roots  of 
the  determinantal  equation  for  a  sequence  of  values  of  the  amount  of  damping 
present.  Here  the  results  of  reference  [6]  will  be  verified  by  applying  the 
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proposed  method  for  analyzing  damped  systems. 

The  equation  of  free  vibration  of  the  system  of  Fig.  T  is 


‘21  O' 

M  - 

i — 

o 

o 

0  21 
L  J 

W 

0  2cJ 

where  I,  c  and  K  stand  for  the  mass  moment  of  inertia  of  the  flywheel,  the  damp 
ing  coefficient,  and  the  torsional  stiffness  of  the  shaft,  respectively,  and 
0i  and  02  represent  the  two  rotational  degrees  of  freedom  of  the  system. 
Equation  (27)  clearly  indicates  that  conditions  (7)  and  (16)  are  not  satisfied. 
Thus,  this  system  does  not  possess  classical  normal  modes  and  the  only  way  to 
determine  its  critical  damping  surfaces  is  through  the  general  method. 

The  determinantal  equation  (4)  for  this  system,  in  view  of  (27),  takes 
the  form 

A“  +  2aA3  +  6AZ  +  6aA  +  5  =  0  ,  (28) 

where 


A  =  A/wx  ,  a  =  c/2Iwt  ,  wj  =  K/I  .  (29) 

For  the  case  of  overdamping  or  critical  damping  for  which  A  =  -B  =  -b/w^  , 

(28)  becomes 

B“  -  2ctB3  +  6B2  -  6aB  +  5  =  0  .  (30) 

Differentiating  (30)  with  respect  to  B  and  requiring  that  (da/dB)  =  0,  one 
obtains  from  (30) 

a  =  2B(3+BZ)/3(1+BZ)  .  (31) 

With  (31),  (30)  takes  the  form 

B6  +  3B*  +  3B2  -  15  =  0.  (32) 

The  only  one  real  positive  root  of  (32)  is  B  =  1.233  for  which  (31)  yields 

the  equation  of  the  only  existing  here  critical  damping  surface  -  a  point,  as 
acr  *  1.474,  (33) 

a  value  which  is  identical  to  that  obtained  in  [6]. 

The  existence  of  just  one  critical  damping  surface  for  this  two  degree  of 


freedom  system  indicates  that  this  actually  corresponds  to  partial  critical 
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damping.  For  0  <  a  <  1.474  the  system  is  underdamped,  for  a  =  1.474  it  is 
partially  critically  damped,  and  for  a  >  1.474  it  is  partially  underdamped. 

Thus,  no  matter  how  great  the  damping  parameter  a  may  be,  there  is  no  way  to 
reach  complete  overdamping,  exactly  because  there  is  no  other  critical  damping 
surface  for  a  >  1.474.  This  observation  was  also  made  in  [6]  by  studying  the 
behavior  of  the  roots  of  (28)  as  the  damping  parameter  a  was  increasing.  In 
fact,  it  was  also  shown  in  [6]  that  for  a  approaching  to  infinity,  (28)  possesses 
only  two  imaginary  conjugate  roots  which  is  physically  explained  by  the  fact 
that  at  a  +  ®  the  damped  flywheel  becomes  locked  and  the  other  flywheel,  with 
its  two  shafts,  forms  a  one  degree  of  freedom  conservative  system. 

This  example  clearly  demonstrates  the  simplicity  of  the  proposed  method 
for  determining  critical  damping.  In  comparison  reference  [6]  required  a  labor¬ 
ious  numerical  computation  and  detailed  plot  of  the  real  and  imaginary  parts  of 
the  four  roots  of  equation  (28)  for  a  lengthy  sequence  of  values  of  the  damping 
parameter  a. 

Example  3 

Consider  the  two  degrees  of  freedom  spring-dashpot-mass  system  shown  in 

Fig.  2,  where  K.  ,  c.  and  m.  (i  =  1,2)  stand  for  spring  constants,  coefficients 

of  viscous  damping  and  masses,  respectively,  and  x  and  x  represent  the  two 

1  2 

translational  degrees  of  freedom  of  the  system.  Notice  that  for  c^=  0  one  has 
a  case  similar  to  that  treated  in  Example  2.  The  equations  of  free  motion  of 
this  system  are 


For  this  system,  equation  (4)  becomes 

X1*  +  2(8  +  8  )X3  +  (w2  +  w2  +  w2  +  48  8  }X2  +  2(8  w2  + 
12  1  2  12  12  12 

+  8  w2  +  8  w2  )  X  +  w2w2  =  0  , 

2  1  2  12  12 


(35) 
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where 


w2  *  K  /m  ,  w2  *  K 
1112  2 

0  *  c  /2m  ,  0 

1  l  2 


111  22 
while  equation  (9)  takes  the  form 


/m  ,  w2  *  K  /m  ,  1 

2  12  21  l 

c  /2m  ,  J 


(36) 


b1*  -  2(0  +0  )b3  +  (w2  +  w2  +  w2  +  40  0  )b2  -  2(0  w2  + 
12  1  2  12  12  12 

+  0  w2  +  0  w2  )  b  +  w2w2  *  0  . 

21  2  12  12 

The  conditions  of  critical  damping  described  by  (10)  yield 

b3  -(3/2) (0  +0  )b2  +  (1/2) (w2  +  w2  +  w2  +  40  0  )b  - 
12  1  2  12  12 

-  (1/2) (0  w2  +  0  w2  +  0  w2  )  *  0  . 

12  21  2  12 


(37) 


(38) 


Thus,  the  critical  damping  surfaces  of  this  system  are  curves  in  the  (0  ,0  ) 

1  2 

plane  described  by  (37)  with  b  =  h.„  being  a  function  of  0  and  0  which 

cr  12 

can  be  obtained  by  solving  (38).  Construction  of  these  curves  can,  in  general, 
be  accomplished  by  computing  the  left  hand  sides  of  (37)  and  (38)  for  all  pos¬ 
sible  combinations  of  a  finite  number  of  values  for  0  ,  0  and  b,  and  select- 

1  2 

ing  those  combinations  that  simultaneously  satisfy  (37)  and  (38);  these  com¬ 
binations  provide  the  points  (3lCr»Szcr)  of  the  critical  damping  curves.  In 
this  case,  however,  the  construction  of  the  critical  damping  curves  can  be 
done  more  easily  by  utilizing  well  known  formulae  (e.g.  [7]  p.  86)  to  solve 

explicitly  the  cubic  algebraic  equation  (38)  for  b  in  terms  of  0  and  0  . 

1  2 

In  general,  equation  (38)  may  have  one  reaT  root  or  three  unequal  real  roots 


depending  on  the  values  of  0  and  0 

1  2 


In  this  problem  one  is  interested  in 
positive  roots  of  (38)  which  upon  substitution  into  (37)  satisfy  it  for  real 


positive  values  of  0^  and  0^  .  One  can  very  easily  prove  by  utilizing  (37)  that 

the  slope  d0  /d0  of  the  curves  at  the  points  0  =0  ,  0  a0  and  0  =0  is  negative. 
2  1  12  12' 

Fig.  3  shows  the  two  critical  damping  curves  C  and  C  of  the  two  degrees 


of  freedom  system  of  Fig.  2  in  the  (0,0)  plane  for  the  following  numerical 

1  2 
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data: 


:  =  3000  lb/ft  =  43,779.528  N/m, 

;  *  4000  lb/ ft  =  58,372.703  N/m, 

2 

i  *  m2  =  1  lb.sec2/ft  =  14.593  N  sec2/m. 


(39) 


(1  ft  =  0.3048  m,  1  lb  =  4.448  N). 

The  two  critical  damping  curves  were  constructed  with  the  aid  of  a  computer 

by  both  the  purely  numerical  method  and  the  one  based  on  the  solutions  of 

the  cubic  equation.  It  was  found  that  for  the  curve  the  cubic  equation 

had  only  one  real  root,  while  for  the  curve  C  it  had  three  real  roots  from 

2 

which  only  one  was  acceptable.  For  8  *  8  =8,  condition  (16)  is  satisfied 

1  2 

with 


a  *  28,  a  =0 
1  2 

and  thus  (21),  (22)  and  (40)  lead  to 

bicr  3  ®cr  ~  wo i  *  ^  ~  ^ 

where  the  w  .'s  are  the  positive  real  roots  of  the  frequency  equation 


(40) 

(41) 


w4  -  (w2  +  w2  +  w2  )  w2  +  w2  w2  =  0  ,  (42) 

O  x  2  120  12 

which  can  be  obtained  from  (35)  for  8  =8  =  0  and  X  =  jw  .  For  the  data  of 

l  2  o 

(39),  equation  (42)  yields 

w  =  24.78  rad/sec  ,  w  =  69.90  rad/sec.  (43) 

Ox  02 

The  satisfaction  of  (16)  also  verifies  the  fact  that  q  =  n  »  2  for  this  system. 

The  critical  damping  curves  C  and  C  divide  the  plane  (8,8)  into  the 

12  12 

regions  R  ,  R  and  R  ,  as  shown  in  Fig.  3.  Curve  C  ,  for  which  8  *8  **1  =  24.78, 
123  1  1  2  Oj 

defines  a  state  of  partial  critical  damping  and  separates  a  region  of  complete 

underdamping  (region  R  )  from  a  region  of  partial  underdamping  (region  R  ). 

1  2 

Curve  C  ,  for  which  8  =8  =w^  =69.90,  defines  a  state  of  complete  critical  damp- 

2  1  2  02 

ing  and  separates  a  region  of  partial  underdamping  (region  R  )  from  a  region 

2 

of  complete  overdamping  (region  R?).  The  results  of  Fig.  3  were  also  verified 
numerically  by  determining  on  the  computer  all  the  four  roots  of  equation  (35) 


with  the  aid  of  Newton's  iterative  method  (e.g.  [5])  for  all  the  possible  com¬ 
binations  of  numerical  values  of  and  6^  ranging  from  0  to  120  in  steps  of 
5.  Indeed,  the  results  indicated  that  all  the  roots  in  region  R^  were  complex 
conjugate  (complete  underdamping),  all  the  roots  in  region  were  negative 
real  (complete  overdamping) ,  while  two  roots  were  negative  real  and  two  complex 

conjugate  in  region  R  (partial  underdamping).  Furthermore,  for  the  point  8  = 

2  1 

8^*6  on  the  curve  there  were  found  two  complex  conjugate  roots  and  two 

negative  real  roots  both  equal  in  magnitude  to  wq^  3  24.78  (partial  critical 

damping),  while  for  the  point  6  *8  *8  on  the  curve  C  there  were  found  two 

12  2 

negative  real  roots  both  equal  in  magnitude  to  w  *  69.90  and  two  more  unequal 

02 

negative  real  roots  (complete  critical  damping). 

The  interesting  thing  in  this  example  is  the  fact  that  region  R^  is  un¬ 
bounded  and  that  even  for  values  of  8t  or  8^  (but  not  both)  approaching  infinity 
there  is  still  partial  underdamping,  while  one  might  expect  to  achieve  over¬ 
damping  if  one  of  the  8's  Is  large  enough.  This  phenomenon  is  the  two-dimensional 
counterpart  of  that  of  Example  2  and  is  amenable  to  an  analogous  physical  ex¬ 
planation.  Indeed,  at  the  limit,  as  8  »  for  example,  approaches  infinity,  the 
mass  m^  becomes  essentially  locked,  x^  approaches  zero,  the  system  becomes  a 
single  degree  of  freedom  one,  and  (34)  reduces  to 

x+28x+w2x=0,  (44) 

2  2  2  2  2 

for  which  the  condition  for  overdamping  or  critical  damping  is 

b2  -  28  b  +  w 2  3  0  .  (45) 

2  2 

Use  of  (10)  then  yields 

8icr  =  w*  *  44'72  •  (46> 

Similarly,  at  the  limit,  as  8  approaches  infinity,  mass  m  becomes  locked, 

2  2 

x  approaches  zero,  the  system  becomes  a  single  degree  of  freedom  one,  and  (34) 

2 

reduces  to 


x  +  28  x  +  (w2  +  w2  )x  *  0  , 

•  I  I  ’  4  1  A  I 


(47) 
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for  which  finally  (10)  yields 


8  „  3  Ar  + »r  =  59.16  .  (48) 

icr  i  u  v 

Equations  (46)  and  (48)  are  equations  for  the  asymptotes  of  the  critical  damping 

curves  and  indicate  that  even  for  infinite  values  of  B  (B  )  there  are  values  of 

1  2 

B  (8  )  in  region  R  ,  which  are  less  than  their  critical  values,  for  which 
2  1  2 

there  is  partial  underdamping. 

During  the  construction  of  the  critical  damping  curves  C  and  C  it  was 

i  * 

observed  that  the  values  of  b  „  and  b  „  were  very  close  to  w  and  w  , 

icr  2cr  J  oio2 

respectively,  for  any  point  (B  ,B  )  on  the  curves.  This  means  that  one  can 

1  2 

approximately  construct  the  curves  C  and  C  by  assuming  that  (41)  holds  true 

1  2 

not  only  for  the  points  8  =8  =8  of  the  curves  but  for  any  point  of  them.  Thus, 

1  2 

(37)  with  b.  =  w  .  finally  becomes 
zk»  r  oz 

A. 8  -  B. 

a  _  1 _ 2-  Ia  o\ 


£i8.  -  h 


where 


A-  3  2w  .  (w2  +  w2 .  )  , 

Z  OZ  2  OZ 

B.  *  w*.  +  w2  w2  +  w2.  (w2  +  w2  +  w2  )  , 

I  Ol  12  Ot.  1  2  12 

£  ■  -  4w  .  <“> 

l  01- 

D.  =  2w  .  (w2  +  w2  +  w2.)  ,  i  -  1,2  . 

X  Ot  '  1  12  ot 

It  is  evident  from  (49)  that  the  approximate  critical  damping  curves  are  two 
equilateral  hyperbolas  with  two  branches  each.  These  branches  are  symmetric  a- 
bout  the  intersection  point  of  their  asymptotes  6  =  D./E.  and  8  *  A./E..  Out 

1  Z  Z  2  Z  Z 

of  the  two  branches  consituting  every  hyperbola  only  one  branch,  that  for 
which  8t  *  Bz  3  wq£,  represents  a  critical  damping  curve.  It  is  easy  to  deter¬ 
mine  the  sign  of  the  slope  of  these  branches  by  computing  dB  /dB  from  (49) 

2  1 

and  utilizing  relations  (50).  Thus,  one  has 

(dB  / dg  )  *  -4w2.w2w2  /(E.B  -  D.)2  , 

21  01212  1  1  1 


(51) 


Indicating  that  the  slope  Is  negative  everywhere. 


A  plot  of  the  approximate  critical  damping  curves  on  the  basis  of  the 
numerical  data  (39)  is  shown  In  Fig.  3.  The  approximate  curves  of  course 
coincide  with  the  exact  curves  at  6  *  6  ,  and  they  slightly  depart  from  them 

i  2 

more  and  more  as  the  B's  increase.  Such  approximate  critical  damping  curves 
can  be  very  useful  for  rapid  design  calculations,  although,  unfortunately, 
no  general  extensions  of  their  validity  can  be  presented  here. 

i 
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Conclusions. 

The  following  remarks  can  be  made  on  the  basis  of  the  preceding  develop¬ 
ments  about  the  free  vibrations  of  linear  discrete  structural  systems  with 
viscous  damping  varying  among  their  elements: 

1)  There  exist  "critical  damping  surfaces"  for  every  system  representing  the 
loci  of  combinations  of  damping  leading  to  partial  or  complete  critically 
damped  motion  and  thus  separating  regions  of  partial  or  complete  under¬ 
damping  from  those  of  overdamping. 

2)  A  general  method  is  proposed  for  determining  the  equations  of  these  critical 
damping  surfaces.  The  determination  of  the  surface  point  corresponding  to 
equal  amounts  of  damping  is  considerably  simplified  for  systems  which,  on 
the  assumption  that  all  amounts  of  damping  are  equal,  possess  a  damping 
matrix  of  the  Rayleigh  type. 

3)  The  dimension  of  these  critical  damping  surfaces  is  equal  to  the  number  of 
Independent  damping  values  in  the  system,  while  their  number  is  less  or 
equal  than  the  number  of  the  system  degrees  of  freedom.  For  a  n  degrees 
of  freedom  system  with  a  damping  matrix  reducible  to  Rayleigh  type  for 
equal  amounts  of  damping,  there  are  at  most  one  completely  underdamped 
region,  one  completely  overdamped  region  and  n-2  partially  underdamped 
regions. 

4)  More  research  is  needed  in  connection  with  the  determination  and  identifi¬ 
cation  of  the  critical  damping  surfaces  for  large  order  systems  which 
require  an  efficient  numerical  treatment. 
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Figure  1:  The  two  degrees  of  freedom  system  of  Example  2. 
Figure  2:  The  two  degrees  of  freedom  system  of  Example  3. 
Figure  3:  Critical  damping  curves  for  the  system  of  Figure  2. 
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